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ABSTRACT 

We examine deep far-ultraviolet (1600 A) imaging of the Hubble Deep Field- 
North (HDFN) and the Hubble Ultra Deep Field (HUDF) to search for leaking 
Lyman continuum radiation from starburst galaxies at z ~ 1.3. There are 21 
(primarily sub-L*) galaxies with spectroscopic redshifts between 1.1 < z < 1.5 
and none are detected in the far-UV. We use two methods to estimate limits 
to the escape fraction (/esc) for these galaxies. First, to compare with previous 
works, we assume a fixed 1500 A to Lyman continuum ratio, /1500//7005 intrinsic 
to the stellar SEDs to convert our /700 limits to relative escape fractions. Second, 
we fit stellar population templates to the galaxies' optical/near-infrared SEDs to 
determine the starburst age and level of dust attenuation for each individual 
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galaxy, allowing a more accurate estimate of the intrinsic /isoo/Ztoo- We find 
that the two techniques do not differ significantly if the galaxies are currently 
undergoing star-formation. We show that previous high-redshift studies may have 
underestimated the amphtude of the Lyman Break, and thus the relative escape 
fraction, by a factor ~ 2. Once the starburst age and intergalactic HI absorption 
are accounted for, 18 galaxies in our sample have limits to the relative escape 
fraction, fesc,rei < 1-0 with some hmits as low as fesc,rei < 0.10 and a stacked 
limit of fesc,rei < 0.08. This demonstrates, for the first time, that most sub-L* 
galaxies at high redshift do not have large escape fractions. When combined 
with a similar study of more luminous galaxies at the same redshift we show 
that, if all star-forming galaxies at z ~ 1 have similar relative escape fractions, 
the value must be less than 0.14 {3a). We also show that less than 20% (3cr) of 
star-forming galaxies at z ~ 1 have relative escape fractions near unity. These 
limits contrast with the large escape fractions found at 2 ~ 3 and suggest that 
the average escape fraction has decreased between z ^ 3 and 2; ~ 1. 

Subject headings: cosmology: observations — galaxies: evolution — ultraviolet: 
galaxies 



1. Introduction 

Both active galactic nuclei (AGN) and young massive stars produce Lyman continuum 
radiation capable of ionizing hydrogen in the surrounding Intergalactic Medium (IGM). The 
fraction of this photoionizing continuum which escapes into the IGM, hereafter referred to 
as the escape fraction (fesc), is a key parameter in determining how both star- formation 
and black hole growth drive the evolution of the ionization state of the IGM. In particular, 
an accurate assesment of the average fesc will help determine if star-formation, rather than 
active galactic nuclei (AGN), is responsible for hydrogen reionization at 2; > 6. 

For the most lumin ous, unobscured AG N (QSOs), the majority of the Lyman continuum 



escapes into the IGM (iTelfer et al.l |2002| ). Since QSO luminosity functions and average 
spectral energy distributions (SEDs) are well known, their contribution to the intergalactic 
photoionizing continuum can be accurately estimated. However, the contribution to the 
ionizing background from star-formation is far more difficult to determine because large 
reservoirs of HI and dust within the galaxy absorb the bulk of the Lyman continuum. 



Studies of the mean transmission through the Lya fores t have measured the total HI- 
ionizing background in the IGM as a function of redshift (IMcDonald fc Miralda-Escude 
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2001: Scott et al. 2002; Bolton et a.. 



2005 



Fan et al.ll2006r) . Recent measurements of the 



QSO luminosity functions at z < 2 ( Richards et al. 20051 ) indicate that QSOs alone provid e 



the vast majority of the intergalactic HI ionizing flux at these redshifts (llnoue et al.ll2006l ). 
At z > 2.5, the QSO space density falls off rapidly and recent measurements of the QSO 
luminosity functions at 2; ~ 3 indic ate that star-forrnation may contribute at least as rnuch t o 
the ionizing background as QSOs (IHunt et al.ll2004j : ISiana et al.ll2007l : iFontanot et al.ll2007l ). 
Since the QSO space density falls off faster than the star-formation density, it is inferred that 
young massive stars mus t provide the ionizing flux responsible for reionizing the universe by 
z ~ 6 Jlnoue et allbooeh . 



The first attempts to measure fp,,r examine d low-redshift galaxies with space based 
ultraviolet (UV) telescopes. iLeitherer et al.l (Il995l ) took UV spec tra of four local star bursts 
which, after accounting for absorpt ion by HI in the Milky Way (IHurwitz et al.l 119971 ) . give 
upper limits of fesc < 0.032 — 0.57. iDeharveng et al.l (1200 ll ) took a deep FUSE spectrum of 



Mrk 45 and derived an upper limit of f^sc < 0.06. 



At higher redshift. iMalkan et al.l (120031 . hereafter M03) took deep HST/STIS broadband 
imaging (A ~ I6OOA) of the Lyman continuum of 11 starbursts at 2 ~ 1.2 and derived 
upper limits to the escape fraction (here, fesc,rei, defined relative to the dust-attenuated 



continuum flux at 1500 A; see Section [2]) fesc,rei ^0.10 — 0.40. ISteidel et al.l ( 200ll . hereafter 
SOI) produced a composite spectrum of 29 Lyman Break Galaxies (LBGs) at z ^ 3.4 and 
detected significant signal below the Lyman continuum. The measured L(1500A) / L(900A) = 
4.6 ±1.0 implies that there is little to no absorption of the Lyman continuum beyond what is 
seen in the UV {lamhdarest ~ 1500A), indicating very small column densities of HI along the 
line-of-sight through the ISM of these LBGs. The authors note that, because of a selection 
bias, this sample is among the bluest quartile of the LBG population and is therefore younger 
and/or less dust-reddened than the average LBG. How ever, other measureme nts at z ~ 3 
using deep optical spectroscopv ( Giallongo et al. 2002[). na rrow-band imaging ( Inoue et al. 



2OO5I ) , and broad-band imaging (iFernandez-Soto et al.ll2003l ) , have all resulted in upper limits 
for fesc,rei between 0.04 and 0.72. 

The first dete c tion o f escaping Lyman continuum from an individual galaxy was reported 
by iBergvall et al. J2006h with a 16ks FUSE spectrum of Haroll, giving an escape fraction 
0.04 < fesc < 0.1, or a relative escape fraction of 0.17 < fesc,rei < 0.42 (though see Grimes 
et al. (2007, submitted) who see no detection down to fesc < 0.02 with a re-reduction of the 
same data). Haroll is a luminous, metal poor, blue compact galaxy - attributes which are 
similar to LBGs at 2; ~ 3. However, it is somewhat underluminous {Luy = 0.2L2bc), may 
have a significantly larger halo mass than typical high-redshift LBGs, and was preselected as 
having a very low HI gas mass M[HI) ~ lO^M© (unpublished result noted in iBergvall et al. 
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20061). 

Most recently, IShapley et al.l (l2006l . hereafter S06) published deep optical spectra of 
an additiona l 14 L BGs with an average fesc,rei ~ 0.14, about 4.5 times lower than the 
Steidel et al.l (120011 ) result. Two of the galaxies appear to have fesc,rei ~ 1, while the remain- 
ing twelve are undetected. This suggests that 10-20% of LBGs have large escape fractions, 
or that the covering fraction of neutral HI and / or dust is 80-90% in these systems. This may 
be the case if mergers or supernovae chimneys (jFujita et al.ll2003l ) sufficiently disturb the HI 
reservoirs, producing low-column density lines-of-sight (LOSs) through the galaxies. 

With only three individual detections of escaping Lyman continuum, it is difficult to 
determine whether any physical parameters (starburst age, lumino sity, dust, HI mass, m or- 
phology) correlate with f^sc- For example, the two detections from IShapley et al.l (120061 ) do 
not stand out from the rest of the sample and are, in fact, very different from one another. 
Furthermore, most programs to date have necessarily targeted the more luminous starbursts, 
so the faint end of the luminosity function has not been sufficiently examined. 



We use deep far-UV imaging of the Hubble D eep Field North (HDF N, IWilliams et al. 



1996|) and the Hubble Ultra Deep Field (HUDF, iBeckwith et al.l 120061 ) to search for the 
escaping ionizing flux from moderate redshift galaxies. The observations are taken with the 
f25QTZ filter on the Space Telescope Imaging Spectrograph (STIS) and the F150LP filter 
on the Solar-Blind Channel (SBC) of the Advanced Camera for Surveys (ACS) aboard the 
Hubble Space Telescope (HST). The F25QTZ and F150LP broadband filters have effective 
wavelengths of ~1600A, providing limits on the Lyman continuum luminosities of 24 star- 
formin g galaxies at z ^ 1.3 . These galaxies are t ypically less luminou s than those of previous 



works (ISteidel et al.l 1200 ll : iMalkan et al.l l2003l : IShapley et al.l |2006|) , allowing us to probe 
trends with luminosity. In Section [2] we examine the various assumptions made about galaxy 
properties which affect our derivation of the escape fraction. Section [3] explains the data set 
and sample selection. In Section 14.11 we use a simple method of applying these assumed 
corrections to the /1500//700 fiux ratio to estimate the escape fraction. 



In Section 14.21 we 

fit stellar population models to the galaxies' SEDs to determine the starburst age and dust 
attenuation. In Section El we compare the differences of these two methods and compare 
our limits to other surveys. 



Throughout the text, we use a flat ACDM Cosmology with Hq = 70 km s ^ 
= 0.3, and = 0.7. All flux ratios use units of flux per unit frequency. 



Mpc" 
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2. Determining the Escape Fraction 

There are two commonly used definitions of the escape fraction in the literature. The 
absolute escape fraction, fesc, is simply the fraction of emitted Lyman continuum photons 
(typically measured at 900A) that escapes into the IGM. This definition is useful in the- 
oretical or semi-analytic models where it is used to directly translate star-formation rates 
and initial mass functions (IMFs) into ionizing backgrounds. Observationally, the absolute 
escape fraction is difficult to determine because the intrinsic Lyman continuum fiux must 
first be known. At low redshift, the Ha line fiux (corrected for reddening by measuring the 
Balmer decrement) is often used to approximate the number of Lyman continuum photons 
emitted from massive stars. However, at high redshift this becomes increasingly difficult as 
the Balmer lines shift to the near-IR. Therefore, the Lyman continuum fiux is often inferred 
from the galaxy's rest-frame UV and optical photometry, in particular the rest-frame 1500A 
fiux, /^(1500A), since that is easily measured aX z > 1. Unfortunately, the 1500A fiux is 
subject to significant dust attenuation so this must be accurately accounted for in order to 
use it for determining the intrinsic Lyman continuum fiux. 



Steidel et al.l (120011 ) defined the relative escape fraction, fesc,rei as the fraction of escap- 
ing Lyman-continuum photons divided by the fraction of escaping photons at 1500A. This 
definition doesn't require knowledge of the level of dust attenuation, and is useful at high 
redshift because /1500 is easily measured and is commonly used when determining luminosity 
functions. Therefore, the relative escape fraction can be used to directly convert luminosity 
functions to ionizing backgrounds. 

In practice, fesc,rei is determined by comparing the Lyman continuum fiux with the 
rest-frame 1500A fiux, /isoo- Spectroscopic and narrow-band studies are sensitive to Lyman 
continuum fiux at wavelengths just short of the Lyman Break (~ 900A), whereas broadband 
photometric observations probe the Lyman continuum at somewhat shorter wavelengths 
(typically ~ 700 A). This observed fiux ratio between 1500A and the Lyman continuum is 
affected by several factors and can be expressed as the following product 



ifimo/hcUs = {fimo/fLc)stei X 10"°-^(^-o"-^^^) xexp(rH/-/GM(i^C^)) xexp(r^,,_,5M(^C)) 

(1) 

where LC is the wavelength at which the Lyman continuum is being observed (~ 700A in 
this study), {fi5oo / fLc)stei is the intrinsic flux ratio from the SED of the stellar population, 
(y4i5oo— ^Lc) is the differential dust attenuation (in magnitudes), thi-igm{LC) is the optical 
depth of the Lyman line and continuum absorption through the IGM along the line-of-sight 
to that galaxy, and thi^isai{LC) is the optical depth of the Lyman continuum absorption 
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from HI within the observed galaxy's interstellar medium (ISM). 

We can rearrange this equation to give us the relative escape fraction 



fescrel = ^^^^^^/^^ X exp{THI-IGM{LC)) = exp{-THj_jsM{LC)) X IQ^-^i^^^^o.-A^a) _ (2) 
1/1500/ jLCjobs 

If we can estimate the amplitude of the intrinsic stellar Lyman break, {fi5oo / fLc)steh and 
the average optical depth of the Lya forest, thi-igm{LC), then then relative escape fraction 
can be computed directly from the observed flux ratio, (/i5oo//lc)o6s- In this section we 
detail each of these factors and quantify their effects on the (/1500 //lc)o6s ratio. Note 
that the relative escape fraction is the product of the HI absorption and the differential 
reddening between 1500A and the Lyman continuum, whereas the absolute escape fraction 
is the product of the HI absorption and the total dust attenuation of the Lyman continuum. 
Therefore, the absolute escape fraction is equal to the relative escape fraction times the dust 
attenuation at 1500 A, 



f _ in-0.4A(1500) f (n\ 
J esc J esc,rel \'^) 



where A(1500) = 10.33E(5 - V) for a Calzetti reddening law (jCalzettil 119971 ). 



2.1. Stellar Population 



The intrinsic flux decrement across the Lyman Break is dependent upon the age, star 
formation history (single burst, exponential decay, constant), initial mass-f unction (IMF), 
and m etallicity o f the stellar population. W e have used both the Starburst99 (ILeitherer et al. 
1999h and BC03 JSruzual fc Charlotlboosh population synthesis models to quantify how these 
properties affect the size of the Lyman Break in star-forming galaxies. We note that the non- 
LTE models whic h better account for ste l lar winds and line blanketing in the atmospheres 
of massive stars (ISchaerer fc Vaccalll998l : ISmith et al.ll2002l ) do not significantly a ffect the 
HI ionizing continuum between 700 < A < 900A (see Figure 8 in lSmith et al.ll2002l ). 



2.1.1. Star Formation History 



Since the lifetimes of the 0-stars (which emit the Lyman continuum flux) are so much 
shorter than the B- and A-stars that dominate the 1500A flux emission, the /1500//700 Aux 
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ratio is highly dependent upon the starburst age and the star-formation history. In Figure [T] 
we plot the expected flux ratio as a function of age for both the BC03 and Starburst99 models 
with single burst and constant star-formation histories. In the single burst systems, the dying 
0-stars are not replenished with new star- formation, causing the /1500//700 ratio to increase 
rapidly after a few Myrs. In the constant star-formation scenario, early in the burst, the 
dying 0-stars are being replenished, while the B- and A-stars accumulate, slowly increasing 
the /1500//700 ratio. Eventually, the B- and A-stars begin to die as well and no longer 
increase in number. Thus, the flux ratio remains constant after ~300 Myr. The flux ratio of 
a starburst with declining star-formation rates evolves between these two extremes. If the 
observed galaxy is undergoing a secondary burst, there may be signific ant remaining f lux at 
1500 A but only if the bursts are separated by less than a few 100 Myr. IShapley et al.l ( 200ll ) 
find that constant star-formation scenarios fit well to the Lyman Break galaxy population 
with a median age, t^/ = 320Myr. We will show in Section 15.11 that we obtain similar star 
formation histories for our sample. Therefore, we expect the typical intrinsic break amplitude 

to be (/l500//900)stei ~ 6 or {fwoo/ f70o)stel ~ 8. 



2.1.2. Initial Mass Function 



The slope of the high-mass end of the stellar initial mass funtion (IMF) determines the 
relative numbe r of 0-stars to B- and A-stars, and therefore impacts the /1500//700 ratio. The 



Salpeter IMF (ISalpeteiill955l ). commonly used in studies of high-redshift galaxy formation, 
has a mass distribution characterized by a = —2.35, where ^(M) oc M". Many studies of 



( Massev 


1998; 


Scalo 


1998; 


Kroupa 


2001) 



of the IMF can be done at high redshift, but initial indications from model fits to the 
UV spectra of a lensed LBG (MS151 2-cB58) indicate th at the high-mass IMF is nearly 



Salpeter and must extend above 50 M© ( iPettini et al.ll2000l ). Given the agreement of detailed 



local studies and initial indications at high-redshift, we use a Salpeter IMF in our analysis. 
It should be noted, however, that several studies suggest that top-heavy IMFs in some 



high-redshift star-formation may offer a possib 



in present-day cluster ellipticals (jFaber et al 



Gibson fc Matteucci 1997; Thomai 



1999 



e explanation of a-element abundances seen 



1992 



Nagashima et al. 



Worthev et all Il992l: iMatteuccil Il994 



2005h . A top-heavy IMF would 



significantly decrease the expected, UV-to-LC ratio expected with a Salpeter IMF. 
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2.1.3. Metallicity 



The metallicities of galaxies such as those in our sample as well as LBGs, range from 
Z ~ 0.2 — 1.0^0. Over thi s range, the arnplitu de of the Lyman break does not change 
significantly (see eg. Fig. 75. iLeitherer et al.lll999l ). In our subsequent analysis we use solar 
metallicity models. 



2.2. IGM absorption 



Along the line of sight (LOS) to any high-redshift galaxy, there are hundreds of inter- 
vening neutral hydrogen clouds whose Lyman line and continuum absorption significantly 
affects the measured far-UV fluxes of the galaxies. The column density and redshift dis- 
tributions of these HI absorbers have been well studied and have been used to model the 



transmission through the IGM as a function of the observed galaxy's redshift (jMadaulll995 



Bershady et al.lll999l ). We make use of t he best available data on the column density and 
redshift distributions of the Lja Forest (IKim et al.l 119971 ) , Lyr nan Limit Systems (LLSs, 
Storrie-Lombardi et al.lll994j ) and Damped Lya systems (DLAs. IStorrie-Lombardi fc Wolfe 
2000l ) to simulate the effects of IGM Lyman line and continuum absorption on the observed 
fluxes of galaxies at 2 ~ 1. We simulate 1000 sight lines to galaxies at redshift intervals 
Az = 0.05 between 1.0 < z < 1.5. For each sight line we randomly place absorbers with col- 
umn densities and redshifts which are consistent with their empirical distributions. We note 
that this does not account for clustering of the absorbers, which may be important, especially 
near the observed galaxy if it is within an overdensity. We then compute the absorption line 
profiles (velocity width (3 = 30 km s~^) and apply continuum absorption for each absorber 
to compute the total transmission through the IGM as a function of wavelength (See Siana 
et al. 2007 for more details). The mean transmission from 1000 LOSs is plotted in Figure [2J 
At the redshifts of interest, the average transmission of UV light changes from 0.57 to 0.37 
between 1.2 < z < 1.5. Therefore, our measurements are more sensitive at lower redshifts 
as less of the Lyman continuum is being absorbed. However the distribution function of the 
IGM transmission is highly non-Gaussian. Figure [3] shows the distribution of IGM transmis- 
sion (as measured through our filter), indicating an almost bimodal distribution. About 20% 
of the simulated lines-of-sight contain large column density absorbers, resulting in very little 
transmission (< 5%). Conversely, about 40% of the LOSs have no intervening large column 
density absorbers and therefore have greater than 60% transmission. Though the mean and 
median are similar (0.51 and 0.59, respectively at z = 1.3), it is clear that the absorption 
correction will change significantly between galaxies. This distribution makes it difficult to 
draw conclusions from small samples if one can't determine the properties of the interven- 
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ing IGM. We will use this IGM transmission distribution in a Monte-Carlo simulation to 
interpret our measurements (see Section [5^ . 



2.3. Dust 



Although it is not necessary to know the level of dust attenuation to determine fesc,rei, 
it is important for converting between the relative and absolute escape fractions and in 
determining whether HI or dust is primarily responsible for the Lyman continuum absorption. 
If the dust extinction continues to increase at A < 912A then the measured ratio of UV-to- 
LC flux will be even larger than the intrinsic Lyman Break of the stellar SEDs. The level 
of increase is defined by the amount of dust attenuation and the reddening law. Of course, 
the extinction curve below the Lyman limit is impossible to measure due to HI absorption, 
but theoretical models of the ISM dust composition and si ze distribution predict extinc tion 
curves that continue to rise to A ~ 700 A (see Figure 14 in IWeingartner fc Drainell200ll ). In 
Figure H] we plot the differential reddening across th e Lyman Break as a function of optical 
color excess, E(B-V), for a Calzetti reddening law (ICalzettil 119971 ) . The reddening law at 
A < 1200 A is extrapolated with a constant slope {dk{X)/dX, where k{X) = E{B — V)/Ax) 
derived at 1100 < A < 1200 A, approxi n iating the increasing extinction predicted below the 
Lyman limit in IWeingartner fc Draind (120011 ). We note that while this extrapoloation is a 
reasonable assumption at 900A (extinction curves in SMC, LMC, MW all continue to rise 
at shorter wavelengths up to the Lyman limit), the extrapolation is somewhat uncertain at 
700 A. The level of reddening in LBGs typically varies from 0.0 < E(B — V) < 0.4 with a 
median value of E{B — V) ~ 0.15 (jPapovich et al.ll200ll : IShapley et al.ll200ll ). Therefore, the 
dust alone can cause the UV-to-LC ratio to vary widely over this range (additional factors 
of 1-8 at 900 A and 1-20 at 700 A) with typical values of the differential dust attenuation of 
~ 2 and ~ 3 at 900 and 700 A, respectively. At color excesses of E{B — V) > 0.5, we can 
not detect the Lyman continuum at 700 A in this study. 



Observations 



We have obtained far-UV imaging of the Hubble Deep Field North and the Hubble 
Ultra Deep Field at 1600 A in three HST General Observer programs (7410, 9478, 10403). 
The HDF-N w'as obs er ved with STI S thro ugh the FUVQTZ filter covering 1.02 arcmin^ 
JCardner et al.lboooh . iTeplitz et all kood \ imaged most of the rest of the HDF-N (3.77 
arcmin^.) with the F150LP filter on the Solar Blind Channel of the ACS. These two instru- 
ment configurations have very similar throughputs (see Figure [5]), with a lower wavelength 
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cutoff at A ~ 1450 A and a decreasing sensitivity out to 2000 A (1840 A for STIS). Tlie effec- 
tive wavelengtlis for tlie STIS and SBC filters are A = I6OOA and A = 1610 A, respectively. 
The HUDF observations (Siana et al. 2007, in prep) use the same ACS/SBC configuration as 
in the HD F-N and cover 7.77 arc min^, nearly the same area as the NICMOS HUDF treasury 
program (IThompson et al.ll2005l ). 

Both the SBC and STIS use Multi- Anode MicroChannel Arrays (MAMAs), which have 
no read noise and are insensitive to cosmic rays. The primary source of noise is dark current, 
which has two components. The first component is a fairly uniform count rate that doesn't 
change with the detectors' temperature. The second component is a temperature-dependent 
"glow" which arises at T > 25 C and is near the center of the SBC MAMA and one corner of 
the STIS MAMA. As the instrument warms up, the dark current from the "glow" increases, 
decreasing the sensitivity of the observations. Because of this, the images taken at the end of 
a series of observations are significantly less sensitive than those at the beginning. Therefore, 
the dark current changes substantially from pointing to pointing, resulting in significantly 
different sensitivities across the field. In addition, there is significant area where adjacent 
frames overlap. Weight maps have been defined which account for the total exposure time 
and dark count per pixel so that accurate sensitivities can be c omputed. The dark s ubtraction 
and weight map production used the same procedures as in iTeplitz et al.l (120061 ). Figure [6] 
shows the cumulative area for which we expect to obtain a 3a detection of a galaxy of a 
given magnitude within a 1" diameter aperture. 



3.1. Galaxy selection 



In order to avoid contamination from photons with Xrest > 912A, we choose galaxies at 
redshifts large enough to shift the Lyman Limit to wavelengths at which our filter/detector 
configurations are no longer sensitive. This occurs at Xobs > 1915A {z > 1.1) for the STIS 
observations and Xobs > 2000A {z > 1.2) for the SBC observations. The MAMA detectors 
are sensitive to wavelengths longer than A > 2000A but their quantum efficiency is lO^'^ 
times lower at 2000 A than at 1500A and decreases rapidly toward redder wavelengths. 

Our sensitivity to Lyman continuum radiation falls off rapidly for galaxies with z > 
1.5 because 1) only the youngest, least reddened galaxies will have significant flux at A < 
6OOA, 2) the rest-frame Lyman continuum of the galaxy will be diminished by continuum 
absorption from an increasing number of neutral HI clouds along the line-of-sight to the 
galaxy. We therefore only choose galaxies with z < 1.5. The spectroscopi c redshifts are 
primarily from the Team Keck Reds hift Survey (TKRS Wirth et al.l 120041 ) in the HDF- 



N and the VLT/F0RS2 spectroscopy (IVanzella et al.l 1200.4 I2OO6I ) in the HUDF, with a few 
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additional redshifts from lCowie et al.l (120041 ) and lLe Fevre et al.l (120041 ) . In addition to these 
redshift requirements, we also ensure that there are no sources with optical spectra with high- 
ionization emission lines (NeV,NeIII) or with strong nuclear point sources indicative of strong 
AGN activity. Finally, we removed the three reddest [B — V > 0.6) galaxies from our sample 
as their SEDs are indicative of older stellar populations with little ongoing star-formation. 
The 8 galaxies in the HDF-N and 13 galaxies in the HUDF that meet these requirements 
are listed in Table [H 

Recent measurements have shown that the total optical throughput of the SBC at red 
wavelengths is much higher than was measured before launch (3.4 x 10~® vs. 9.5 x 10~^ at 
A = 3500 A, STScI Analysis Newsletter, Nov. 30, 200^. The total system throughput at 
3500 A is still four orders of magnitude less than at 1500 A, but we must be careful that 
the 'leaking' optical light does not significantly affect our sensitive far-UV measurements. 
As a test, we examined the effects of the red leak with two types of galaxies: a 'typical' 
galaxy (igf = 300 Myr, E{B — V) = 0.2) and a 'worst-case' (ie. extremely red) galaxy 
{tsf = 1 Gyr, E{B — V) = 0.5). We shifted the galaxy templates to 2; = 1.2, applied average 
IGM absorption, and determined the ratio of leaking flux density (with Xrest > 912 A) to the 
measured flux density at Xrest > 1500 A. These ratios were 0.003 and 0.002 for the worst-case 
and typical galaxies, respectively. Therefore, we can be confident that our results will not be 
significantly affected if our observed flux ratios are significantly above /700//1500 ~ 0.003. We 
show in Section HIT] that our most sensitive limits are /700//1500 ~0.01, and should therefore 
be unaffected by this leaking optical light. 



3.2. Photometry 



For t he optical and near-IR photometry, we us ed the NIGMOS-selecte d photometric 
catalogs of bickinson et all koO(h in the HDF-N and [Thompson et all J2005h in the HUDF. 
The HUDF near -IR photometry was corrected by AJ = 0.30 and AH = 0.18 as discussed in 
Coe et al.l (120061 ). For both of these catalogs, source detection and isophotal apertures were 
determined from the summed J (FllOW) and H (F160W) band NICMOS images. These 
same apertures were used to extract fluxes from the optical images as well. Each isophote 
was inspected by eye to ensure there was no contamination from nearby galaxies and that 
no galaxy was improperly fragmented into several pieces. 

The UV photometry was extracted using isophotes defined by the optical B-band [Xrest ~ 
I900A), which should directly trace new star-formation at these redshifts. Here we make 



http:/ /www. stsci.edu/hst/acs/documents/newsletters 
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the reasonable assumption that significant far-UV signal will not be detected where there is 
no measurable flux in the deeper B-band images. Because the B-band images are so sen- 
sitive, their faintest isophotes have large areas which contain very little flux. We find that 
decreasing the aperture to exclude the faintest 20% of the galaxy reduces the isophotal areas 
by factors of ~2-3. Therefore, in order to maximize the signal-to- noise of the photometry, 
we use as our extraction aperture the isophote which contains only 80% of the B-band flux, 
increasing our sensitivity by ~ 0.4-0.6 mags. No correction for foreground extinction is 
necessary because both the HDF-N and HUDF have clear lines of sight through the galaxy 
{E{B - y) - O.Ol. lSchlegel et al.lll998h . 



As discussed by ( iGardner et al.l l2000l : iBrown et al.l l2000l : iTeplitz et al.l l2006l ) , the dark 
current is the principle source of noice in these observations. The individual frames were 
weighted by the square of the exposure time and divided by the total dark (primary + 
"glow"). As the dark count scales with exposure time, these weight maps scale linearly with 
the ratio of exposure time to dark rate. Thus, the final weight maps are the square of the 
signal-to-noise ratio for objects fainter than the background. For non-detections, we derive 
3a upper limits from these weight maps. 



4. Results 

We do not detect any of the 21 galaxies in our sample (with S/N > 3) and the 3a 
upper limits to the far-UV flux are given in Table [H The distribution of measured SNR 
(See Figure [7]) is centered around zero as expected. There is one galaxy which has a large 
negative flux, but that is due to poor background subtraction due to edge effects. The other 
sources are not affected by this. For the bluer, more luminous galaxies in our sample, we 
expect to derive significant limits to the escape fraction of ionizing photons. In Section 
14.11 we derive fesc,rei limits by reproducing the techniques of previous high-redshift escape 
fraction studies. Specifically, we calculate flux ratios on either side of the Lyman limit and 
compare to the average intrinsic values to deduce a relative escape fraction. In Section 14.21 
we fit SED models to ascertain dust and age parameters for each galaxy to better determine 
the intrinsic Lyman continuum flux and thus the relative escape fraction. 



4.1. Flux Ratios 



To replicate the same methods used in previous high-redshift escape fraction studies, 
the 3a limits to the (/i5oo//70o)o6s flux ratio were computed and are listed in Tabled! The 



- 13 - 



/i500 value is computed by interpolating the optical photometry and the /700 limit is derived 
from the weight maps which account for total exposure time and dark current per pixel. The 
flux ratios are converted to fesc,rei with Eqn[21 We choose to use {fi5oo/ f 700) stei = 8 for our 
intrinsic stellar flux ratio since the ratio varies from 6-10 between 0.01 < tage < 0.3 Gyr. 
We use the redshift dependent average IGM absorption correction exp{—TjGM{z)) from our 
simulations defined in Section 12.21 

The distribution of fesc,rei limits are plotted in Figure [H] (top panel). Nineteen of the 
galaxies have fesc,rei < 1-0, indicating either increased dust attenuation or HI absorption of 
the Lyman continuum. 



4.2. SED fitting 

In addition to deriving flux ratios, we have fit SED models to the high S/N optical/near- 
IR photometry. In this way, we derive better estimates of the starburst age and dust red- 
dening for each galaxy, rather than assuming average values for the entire sample. We used 
all available optical/near-IR photometry from the deep HST images for the SED fits. These 
data, in addition to measuring the slope of the UV continuum, span the 4000 A break, al- 
lowing us to break the degeneracy of age and dust effects on the UV slope. The observed 
far-UV (1600 A) limits are not used in the fits of the SEDs. 

We fit a suite of nine Bruzual-Charlot models with different starburst ages varying 
logarithmical ly between 0.001 < tage < 10 Gyr and allow dust attenuation (Calzetti Law, 



Calzettilll997l ) to vary as a free parameter. The metallicity was fixed at the solar value and the 



IMF is assumed to be Salpeter (ISalpeteiill955l ). We attempt to fit models with instantaneous 



starbursts, exponentially declining star formation with an e-folding time r = 100 Myr and 
constant star-formation models. 

The results of the fits are given in Table [T] and a few examples are plotted in Figure 
M We find that the best fit model is a constant star-formation model for 90% (19 of 21) 
of the galaxies, confirming that these galaxies, which are bright and blue in the rest-frame 
ultraviolet, are actively undergoing star-formation. The fits to the two other galaxies have 
values that are not significantly better than those of the constant star-formation model. 
Therefore, we have chosen to list only the contant star- formation fits in Tableland use only 
these models in the subsequent analysis. 

The distribution of relative escape fractions using the {f 1500 / f 700) stei from the best-fit 
model is shown in Figure [8] (bottom panel). Eighteen of the galaxies have a relative escape 
fraction less than unity. We stacked the fluxes at 1500 A and added the 700 A errors in 
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quadrature to achieve a summed 3a limit fesc,rei < 0.08 at 700 A, though most of the signal 
comes from our brightest source. 



As a consistency check, we also fit to the Starburst99 (ILeitherer et al.lll999l ) models with 
constant star-formation, solar metallicity, and Salpeter IMF. There were small differences in 
age and dust determinations, but no systematic difference from the results derived with the 
Bruzual & Chariot models. 



4.3. Individual Galaxies 

4.3.1. J123652.69+621355.3 

This is the brightest source {B = 22.41) in our sample and is undetected in the far-UV. 
However, there is a clear {7a) detection of a compact source ~ 0.3" to the North (Figure 
[TOl) . This smaller object is part of the larger isophote used in computing the optical/near-IR 
fluxes but is not in our smaller aperture used for far-UV photometry. It has a flat spectrum 
(in f^) across the UV/optical and a far-UV flux /700 = 0.030 /iJy, with no break in the SED 
{FUV — U (AB)= 0.07). Therefore, it is unlikely to be at the same redshift since there is no 
indication of IGM absorption. However, it is interesting to note that, if this system was at 
2 ~ 3 and was observed with ground-based (ie. low spatial resolution) optical spectroscopy, 
this object would appear to have escaping Lyman continuum, with fesc,rei ~ 0.20. Though 
this is the only object in our sample which appears to have foreground contamination, the 
contamination at higher redshifts will be much higher as the foreground path is larger. 
Therefore, it may be important to obtain high-resolution follow-up of high-redshift escape 
fraction detections to ensure that the flux is not originating from low-luminosity foreground 
objects. 



4.3.2. J033239. 09-274601. 8 



There is a bright {B = 21.4) QSO at z = 1.22 fKX4. ICroom et all 120011 ) which was 
not included in this analysis but is detected at ia in the far-UV with /700 = 0.015 ± 0.004 



/iJy. The optical/near-IR photometry is best flt with power law slope a 



-0.85 



where 



fu oc z/""), consistent with the average near-UV spectral slope of QSOs ( iTelfer et al. 



2OO2I ). 



If we extrapolate this power-law and apply a correction for IGM absorption, we expect to 
measure a far-UV flux /700 ~ 0.95 yuJy. Therefore, we estimate an absolute escape fraction 
fesc ~ 0.02, a very small value for an unobscurred QSO. It is also possible, however, that the 
escape fraction is high, but that a high column density absorber lies along the line of sight 
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to the QSO, as is expected for ~ 20% of the LOSs at this redshift (see Figure [3]). 



Discussion 



5.1. Comparison of SED Fits to Flux Ratios 

The best fit ages and color excesses span the allowed ranges, wi th median values of 
t = 300 Myr and E(B -V) = 0.19, similar to the values found in LBGs Jshaplev et al.lboOll : 
Papovich et al.ll200ll ). We used the best fit SEDs to compute the (/i5oo//70o)sie« flux ratio 



and compare with the assumed value of 8 used in section 14. 1[ The intrinsic flux ratio spans 
the range of 6.1-11.2 with a median of 9.6, 20% higher than our assumed average value but 
within the range of uncertainty of other parameters. We conclude that simply correcting the 
measured flux ratio with an average (/isoo/ /too )stei ~ 8 — 10 is not unreasonable for galaxies 
with ongoing star-formation (eg. LBGs). The distribution of relative escape fractions using 
both methods (constant UV-to-LC ratio and UV-to-LC ratio from the SED fits) are plotted 
in Figure [8], and show very similar distributions. In the subsequent analysis of individual 
objects, we use the values derived from the SED fits. 

We have computed the UV-to-LC ratio at 900A, (/i5oo//90o)stei, from our fits to compare 
with the spectroscopic studies, finding a med ian value of 7 . 3. Th is is more than a factor o f 
two larger than the value of ~ 3 assumed by ISteidel et al.l (l200ll ) and IShapley et al.l (120061 ) 
and would imply an a verage relative escap e fraction ~ 2 times higher than those estimates. 
This correction to the Steidel et al. J200lh measurement would yield an unphysical relative 
escape fraction greater than unity. Their sample is from the bluest quartile of the total LEG 
population and may be amongst the youngest, thus decreasing the expected intrinsic break 
(See Figured]). However, it is difficult to imagine how the incidence of star- formation can 
be synchronized across entire galaxies (seve ral kpcs) such that t he the total star-formation 
in the galaxy is less than 10 Myr old. The IShapley et al.l (120061 ) sample is more indicative 
of the total LEG population, where larger values for the intrinsic Lyman break should be 
used when computing the average fesc,rei- However, our fits show a factor ~ 2 variation in 
the break amplitude, demonstrating the need to deduce the starburst age when analyzing 
individual galaxies. 



5.2. HI Column Density Limits 



A fundamental parameter in escape fraction studies is the relative importance of dust 
and HI in attenuating the Lyman continuum. If we assume a Calzetti extinction curve and 
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extrapolate to from 1200A to 700 A where the extinciton curve is not empirically determined, 
we can use our best fit values for dust extinction to determine the level of attenuation at 
700 A. If our measured flux limit is still lower than the expected flux after dust attenuation 
then there must be additional absorption by neutral Hyrdrogen. We can then derive lower 
limits to the HI column density within the ISM of these galaxies. These HI ISM transmission 
limits are given as exp{—THi,isM) in Tabled! Our large lower limits to the observed /1500//700 
ratio (after correcting for the intrinsic break and IGM opacity) can be explained entirely by 
differential dust extinction in all but four of the galaxies. The limits to the HI transmission 
for these four galaxies are 0.45, 0.61, 0.78, 0.88 and are not strong limits to the HI column 
densities. Significantly deeper observations are needed to determine whether dust or HI is 
the principal cause of Lyman continuum absorption within the galaxies' ISM. 



5.3. Foreground Contamination 

One of our galaxies (J123652. 69+621355. 3) has a far-UV detected object within 0.3" 
of the aperture. These two objects would appear as one if observed with lower spatial 
resolution. Therefore, it is possible to have a foreground object contaminate the photometry 
(or spectroscopy) of a high redshift galaxy so that it appears to be emitting in the Lyman 
continuum. This is especially true with ground-based studies where the spatial resolution 
is low. The only detections of escaping Lyman continuum have been found at 2; ~ 3 with 
ground-based studies so care must be taken to ensure that foreground contamination is not 
a problem. 

Galaxies at 2; ~ 1 have larger angular diameters than galaxies at 2 ~ 3, which increases 
the likelihood of contamination. However, the comoving line-of-sight distance to 2; ~ 3 is 
more than twice the line-of-sight distance to z ~ 1 and the space de nsity of UV-luminou s 



galaxies is larger between 1 < z < 3 than it is at lower redshift (lArnouts et al.l l2005l ). 
Therefore some contamination in z ~ 3 studies might be expected. 

We can use galaxy number counts in the f/-band to determine the likelihood of fore- 
ground contamination of 2; ~ 3 as these galaxies must reside in the foreground, < 2 < 3, and 
are emitting at wavelengths that mimic Lyman continuum of LBGs. Because the expected 
(/i50o//90o)obs is SO high, even very faint foreground sources {U ~ 28) can cause L* LBGs to 
appear to have large escape fractions. The surface density of objects with U{AB) < 28 is 
~ 3.5 X 10^ deg~^ (Dolch et al 2007, in prep). If we assume that ground-based imaging or 
spectroscopy can not resolve objects which lie within a 0.5" radius of each other, than we 
would expect that each 2; ~ 3 galaxy has a ~ 2% chance of foreground contam ination. Given 



this probability, there is a ~ 22% chance that one galaxy (out of 14) in the IShapley et al. 
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( 120061 ) study is subjected to foreground contamination, but only a ~ 3% chance that the 
both detections are contaminated. Therefore, it is unhkely that foreground contamination 
can entirely explain the large escape fractions at 2; ~ 3. However, it is important to keep in 
mind that only a small percentage (maybe ~ 10 — 30%) of LBGs appear to exhibit large es- 
cape fractions. Although foreground contamination will only cause a few percent of all LBGs 
to appear to have escaping Lyman continuum, they may comprise a significant percentage 
of LBGs chosen to have high escape fractions (~ 10 — 20%). This demonstrates the need for 
high resolution follow-up imaging of sources to confirm that there is no contamination. 



5.4. Combined Analysis with Other z ^ 1 Studies 



Malkan et al.l (l2003l . herafter M03) conducted a similar study to ours, observing 11 
luminous starbursts at 1.1 < z < 1.5 with STIS and obtained no detections. They use 
a significantly broader and bluer filter (F25SRF2) with a Ag// = 1453A, ~ 150 A shorter 
than our observations. Accounting for the mean redshifts of the surveys, this corresponds 
to Xrest ~ 660 A for M03 and Xrest ~ 715 A for our study. Therefore, the flux ratios of M03 
are subject to larger effects from dust attenuation and star formation history than those 
measured here. Although the STIS camera is less sensitive than the ACS/SBC, the broader 
filter increases the sensitivity so that the M03 /700 limits are close to our derived limits. In 
addition, they targeted more luminous starbursts, resulting in limits to the fesc,rei that are 
~ 2 times lower than ours. The limits in Table 2 of M03 are la limits and are an additional 
factor of two too low, as the radius of the aperture used for the limits was smaller, by a 
factor of two, than the intended 0.5" aperture radius (M. Malkan, private communication). 
After applying corrections for IGM absorption (which M03 did not do), we derive limits of 
0.10 < fescM^^) < 0-42. 

We have fit the M03 galaxies to SED models and verified that the constant star- 
formation assumption of (/i5oo//70o)ste« is appropriate for their sample. The STIS field- 
of-view is small so there are no objects (target galaxy or otherwise) detected in 9 of 11 
pointings. The far-UV images can not be astrometrically aligned with optical data, leaving 
their measurements subject to the intrinsic pointing uncertainties of HST. Because these 
pointing uncertainties are larger than their aperture radii (1" vs. 0.5"), the stacked far-UV 
images of the 11 galaxies spread the light over larger areas and do not give significantly 
better limits than the individual images. 

Combining with the M03 sample, we now have 28 galaxies with fesc,rei ^ 1-0 and 20 
galaxies with fesc,rei < 0.5 at z ~ 1. The M03 sample have UV luminosisites similar to L* 
LBGs at 2 ~ 3, whereas our sample is somewhat fainter (0.1-1.0 L^^fj). These galaxies span 
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nearly two orders of magnitude in luminosity and have a broad range in morphologies and 
starburst ages, yet we see no evidence for large escape fractions at this redshift. 

As shown in Section 12.2^ the transmission through the IGM can vary substantially 
along different lines-of-sight. This will affect our results such that ~ 20% of our limits are 
effectively meaningless as they are looking at opaque lines of sight. However, ~ 60% of our 
limits will be substantially better than the limits given in Table [T] as these lines-of-sight are 
more transparent than average. To better evaluate the effects of this distribution on our 
survey, we've performed a Monte-Carlo simulation where we observe 32 galaxies (21 from 
this sample and 11 from M03) with the same distribution of (/1500/ /too )ofes limits as in our 
surveys. The IGM transmission is chosen randomly from the distribution plotted in Figure 
m We assume an intrinsic break of (/1500//700) ste« = 8 for all galaxies. We then need to 
assume an escape fraction for each galaxy. For this we assume that a fraction Y of star- 
forming galaxies at our redshifts have an fesc,rei = X, and all other galaxies have effectively 
no Lyman continuum transmission. This analysis assumes that this parameter space (X,Y) 
is the same for galaxies with different luminosities. We allow X and Y to vary between 0.0 
and 1.0 and run 10,000 iterations of our observations for each parameter set (X,Y). We then 
determine, for each parameter set, what percent of the time we achieve a null result for all 
32 galaxies. The results are plotted in Figure [TTl The shaded regions denote the parameter 
space which is excluded by our combined samples at 68, 95, and 99% confidence (from lighter 
to darker). If all star-forming galaxies at 2; ~ 1.3 have the same relative escape fraction. 
Figure [TT] shows that it must be less than fesc,rei < 0.14 at 99% confidence. Conversely, if 
some galaxies have large relative escape fractions, fesc,rei > 0.75, while the others have none, 
they must be less than 20% of the total population (at 99% confidence). 



5.5. Uniform Comparison of Previous Studies 



In Table m we summarize the escape fraction studies at all redshifts by converting the re- 
sults to common definitions: the UV-to-LC ratio corrected for IGM absorption, (/i5oo//90o)corr, 
and the relative escape fraction. 

Because observations at Xobs ~ lO OOA are so difficult, the local sample consists of only 
six objects. Once the Milkv Wav HI (jHurwitz et al.l 119971 ) and foreground dust extinction 
(see corrections in iDeharveng et al.ll200ll ) are properly accounted for, and the measurements 



are converted to relative escape fractions, the data do not p ut strong limits o n the relative 
escape fraction. In addition, three of the four galaxies in the iLeitherer et al.l (119951 ) sample 
have very significant color excess, E{B — V) > 0.5. If the Lyman continuum is subject to 
the same dust extinction, then it would not be detectable in these measurements. The most 



- 19 - 



sensitive limits at low redsliift come from FUSE spectra of Mrk 54 flDeharveng et al.ll200ll ) 
and Haro 11 (Grimes et al. 2007, submitt ed) which give UV-to-LC limits of 112 and 21 
respectively. Note that iBergvall et al.l (120061 ) originally reported a strong Lyman continuum 
detection, but a recent reanalysis by Grimes et al. (2007, submitted) detect no flux blueward 
of the Lyman limit. 

The 2; ~ data set is too small to rule out large escape fractions in a small subset of 
galaxies, and the limits are too weak to rule out small escape fractions in a high percentage 
of galaxies. A larger, deeper s tudy of low redshift starbursts (perhaps the UV- luminous 
galaxies of iHoopes et al.l (120061 ) ) is required for better comparison with the high redshift 
studies. 

The z ~ 3 studies claim markedly different conclusions and warrant a more detailed 
discussion. In Table [2], we have converted the results of the studies to common quantities. 



Steidel et al.l (I2OOII ) found a significant (nearly 5a) detection of the Lyman continuum 
at Xrest = 880 — 910 A in a composite spectrum of 29 2; ~ 3.4 LBGs. Once correcting for 
IGM absorption, the implied UV-to-LC ratio, (/i5oo//90o)corr = 4.5 implies a relative escape 
fraction near un ity. Note tha t the s ensitivity could not place significant constraints on any 
individual LBG. ISteidel et al.l (I2OOII ) point out that, because the LBGs were required to be at 
the high redshift end of the ^7-dropout sample, these 29 LBGs were among the bluest quartile 
of the entire LBG sample. Therefore, it might be expected that differential reddening is not 
a major issue, and the intrinsic Lyn i an bre ak amplitude is somewhat smaller (due to younger 
ages). Regardless, as ISteidel et al.l (I2OOII ) point out, this is a surprising result which must 
be verified with deeper spectra. 



Shapley et al.l (120061 ) published much deeper individual spectra of 14 LBGs, of which 



two showed significant flux below the Lyman limit. These UV-to-LC ratios (2.9, 4.5) also 
imply relative escape fractions near unity, while the spectra of the remaining 12 give a 3a 
limits from fesc,rei < 0.25 — 1.0, depending on depth (when using a Lyman break amplitude 
of 6). These limits still allow for non-zero escape fractions, but do not display the small 
UV-to-LC ratios exhibited in the two detection or in the stack of ISteidel et al.l (1200 ll ). 



Several other groups have published escape fraction limits of LBGs at 2; ~ 3, finding 
no detections. Once converting the limits to 3a and using a mo r e con serva tive Lyman 
break amplitude, (/i5oo//90o)stei = 6, the data from iGiallongo et al.l (120021 ) and llnoue et al. 



( I2OO5I ) do not significantly constrain th e relative escape fraction (ie. limits of fesc,rei ~ 1, See 
Tabled). iFernandez-Soto et al.l (120031 ) have analyzed broadband photometry of 27 galaxies 
between 1.9 < z < 3.5 in the Hubble Deep Field-North and report a 3a limit of fesc < 0.039. 
This low limit, derived from such a large sample, appears to contradict the large escape 
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fractions detected by SOI and S06. How ever, there are several differ ences in the analysis 
which should be accounted for. Firstly, the lFernandez-Soto et al.l (120031 ) SEDs do not include 
any intrinsic Lyman Break in the stellar SEDs, so the limits need to be multiplied by the 
additional {fi5oo / hoo) stei = 8 (the value listed in Table [2] accounts for this). Secondly, the 
filter used to measure the escaping Lyman continuum contains significant flux from redward 
of the Lyman break, necessitating very accurate determinations of the UV continuum flux 
in order to estimate the level of Lyman continuum in the same filter. Once these effects 
are accounted for, th eir stacked limit is above the average detection of S06. In summary, 
Shapley et al.l (120061 ) is the only large study at 2; ~ 3 which significantly constrains the 
escape fract ion, and it confirms the preliminary detection of significant ionizing emissivity 
at z ~ 3 bv lsteidel et~aD J200lh . 



The best limits on the relative escape fr action at z ^ 1 (Section 15.41) are signifi cantly 
lower than the average detections at z ^ 3 (ISteidel et al.ll200ll : IShapley et al.ll2006l ). and 
yet we find no detections in the present survey. In Figure [11], we show the parameter space 
implied by the S06 results (dashed box) does not agree, at la, with our results, implying 
that there may be an evolution in the relative escape fraction with redshift. However, many 
of the galaxies in our combined sample are somewhat less luminous than the ^ ~ 3 studies 
{Luv = 0.1 — l.OL^y for z ~ 3 LBGs). It is possible that higher UV lu minosities, and 
thus star-formation rates, can create supernova chimneys (jFujita et al. 2003) which clear the 
surrounding ISM or that the high Lyman continuum flux can ionize low HI column density 
sight lines. Therefore, the fact that our galaxies are less luminous may help explain why 
we see no examples of unity escape fraction. The galaxies in the M03 sample have similar 
luminosities to S06 so one might expect to detect a couple of galaxies with large relative 
escape fractions. However, the M03 sample alone is too small to rule out the parameter 
space of S06. 

Another difference between the 2; ~ 1 and 2 ~ 3 studies is that our broadband obser- 
vations sample shorter wavelengths of the Lyman continuum so they are far more sensitive 
to the star formation history. Our assumption of constant star formation is a best-case 
scenario since this minimizes the amplitude of the {fi5oo/ f 700) stei ratio. If there has been 
any significant decrease in star-formation within t < 10 Myr, /700 will be significantly lower 
than /900, while having little significant effect on the broadband SED redward of the Lyman 
Break. Therefore, it is possible that some of our galaxies will have a larger (/i5oo//70o)stei, 
thus weakening our limits on fesc,rei- 
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6. Summary 

We have examined deep far-ultraviolet (1600A) imaging of the HDF-North and HUDF 
to search for escaping Lyman continuum flux from 21 star-forming galaxies at 1.1 < z < 1.5, 
finding no detections. We account for all factors needed to properly convert our limits to 
relative escape fractions, including intervening IGM absorption and the amplitude of the 
intrinisic Lyman break of the stellar population. 

By fitting SED models to the optical/near-IR photometry, we estimate the star-formation 
age and dust reddening parameters for each individual galaxy. We find that, although there 
is ~ 50% scatter in the intrinsic stellar Lyman Break due to the starburst age, it is reason- 
able to assume a constant {f 1500 / froo) stei ~ 8 or {f 1500 / f 900) stei ~ 6 for galaxies with SEDs 
similar to Lyman Break Galaxies. This value is two times larger than the [finoo/ f 900) stei = 3 
assumed in many previous studies. Assuming a reasonable extrapolation of the extinction 
curve below the Lyman limit, we show that the observed flux decrement at 700 A can be 
attributed to dust attenuation and does not require large column densities of HI within the 
ISM of the galaxies. Deeper observations are required to determine the relative importance 
of dust and HI to the escape fraction. 

We obtain 3a limits better than fesc,rei < 1-0 in 18 galaxies and fesc,rei < 0.5 in 9 
galaxies, with some limits down to fesc,rei < 0.10. Our stacked fluxes give a combined limit 



of fesc,rei < 0.08, similar to the sensitivity achieved by iMalkan et al.l (120031 . corrected with 
common assumptions with this study) for more luminous starbursts at the same redshift. 
This is the first study to achieve these sensitivities on high redshift starbursts which are less 
luminous than typical L*j^bg- These stacks give the deepest escape fraction limits achieved at 
any redshift and demons trate a paucity of ion izing emissivity in most starbursts at z ^ 1.3. 



Including the sample of iMalkan et al.l (120031 ). we now have 28 galaxies at z ~ 1.3 with 



fesc,rei < 1 and 20 galaxies with fesc,rei < 0.5 and no detections. When properly accounting 
for the broad distribution of IGM transmission with Monte Carlo simulations, we conclude 
that, at 99% confidence, less than 20% of star-forming galaxies at 2; ~ 1 have relative escape 
fractions near unity. Conversely, if all of the galaxies have similar relative escape fractions 
then, at 99% confidence, they must have fesc,rei < 0.14. 



When all of the existing z ^ 3 studies are uniformly compared, only the IShapley et al. 



( 120061 ) study gives significant limits to the relative escape fraction and finds large ionizing 
emissivity in at least some (~ 10 — 30%) LBGs. The z ~ 1 studies are sensitive to lower 
escape fractions than those at z ~ 3 (even those with detections), and yet no galaxy at 2; ~ 1 
is seen emitting significant ionizing flux. This disagreement is marginally significant (Icr), 
suggesting a possible decrease in the escape fraction with redshift. We cannot yet rule out 
the possibility that the discrepancy is the result of differeint observing methods which probe 
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different regions of the Lyman continuum (700 A vs. 900 A). 

Further investigation of the evolution of the escape fraction requires larger samples of 
z ~ 1 galaxies that are better analogs of the z ~ 3 LBGs. We have an ongoing program of 
HST FUV imaging and spectroscopy to obtain better limits on the escape fraction at 2; ~ 1 
in the GOODS and COSMOS fields. Other projects, such as stacking of GALEX sources, 
will also greatly improve the current limits. 
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Table 1. Far-UV limits and best-fit SED parameters. 
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Table 2. Compilation of flux ratio limits (corrected for average IGM absorption) and 
relative escape fractions. The IGM absorption is redsliift dependent so, in order to 
facilitate comparison between surveys at different redsliifts, we have accounted for this by 
multiplying the flux ratios by the average IGM transmission at the corresponding redshift. 
Limits are in parentheses and have been convert ed to 3cr. The z 'y res ults also account 
for foreground extinction using values derived in lDeharveng et al.l (120011 ) . A conversion of 

/goo = 1-333 X /700 was used based on the constant star-formation SEDs. A 
{fi50o/ f 900) stei = 6 or (f^^m/f7 m)..M. = 8 ratio has be e n assu med to convert (/i5oo//90o)o6s or 
(/i50o//70o)o6s to fesc,rei- The iFemaudez-Soto et al.l (120031 ) limits have been multiplied by 
eight to account for the intrinsic Lyman break of the stellar population. 
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Fig. 1. — The /1500//700 and /1500//900 intrinsic flux ratios as a function of time since onset 
of star-formation using tlie BC03 (solid) and Starburst99 (dashed) models . 
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Fig. 2. — The average transmission through the IGM of UV hght from galaxies at redshift 1.2, 
1.3, 1.4, and 1.5 (left to right). Each curve is computed from simulations along 1000 lines of 
sight through the Lja forest. The corresponding Lyman limits arc denoted by solid vertical 
lines. As a reference, the transmission curve is plotted for the ACS/SBC F150LP filter used 
in our observations (dot-dashed line). The vertical dotted line is the pivot wavelength of the 
filter curve. 
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Fig. 3. — The distribution of IGM transmission within the F150LP filter for 1000 fines of 
sigfit towards galaxies at z = 1.3. Tfie spike at zero transmission is due to Lyman Limit 
Systems (LLSs) and Damped Lya systems (DLAs) witfiin Az ~ 0.4 from tfie target galaxy. 
Tfie mean and median are 0.51 and 0.59, respectively. 
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Fig. 4. — The differential reddening from 1500 A and the Lyman continuum (defined as 
]^Q-o.4(A(i50o)-A(LC)) ^]^gj.g LQ ig either 900 or 700 A) as a function of dust reddening 
for a Calzetti reddening law extrapo lated to A = 700 A. The dotted vertical line shows the 
median value oi E[B—V) = 0.155 for lShapley et al.l (l200ll ) for Lyman Break Galaxies derived 
from broadband photometry. The tail of the reddening distrubution for LBGs extends out 
to E{B -V) = 0.4. The 700 A flux is affected far more than at 900 A showing that our 
measurements are more sensitive to dust than those measurements just blueward of the 
Lyman limit. 
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Fig. 5.— Total system throughput for the ACS+SBC+F150LP (sohd hue) and STIS+FUV- 
MAMA+F25QTZ (dashed) configurations. The ACS/SBC is nearly 3x more sensitive than 
the STIS configuration. The dashed vertical lines denote the location of the Lyman Limit 
at z=l.l and z=1.2, the low redshift cutoffs for our sample selection for each configuration. 
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Fig. 6. — Detectable area as a function of limiting magnitude (AB, 3a) within a 1" diameter 
circular aperture. The solid line is the total area from all three surveys. The dashed, dotted, 
and dash-dotted lines are the areas of the HUDF, HDFN-STIS, and HDFN-SBC surveys 
respectively. The HUDF survey is much larger, allowing us to detect more objects, while 
the HDF-SBC survey is deeper. We are are sensitive to fainter galaxies if their extraction 
isophotes are smaller. 
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Fig. 7. — Histogram of measured signal-to-noise ratios for the 29 galaxies in our sample. 
There are no detections above SNR > 3. The dashed vertical line is the average < SNR >— 
—0.40. The standard deviation of this distribution, a — 1.7, is slightly larger than expected. 
This is due to small errors in "background" subtraction of a few sources (namely the object 
at SNR — —5) since the dark current is non-planar. 
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Fig. 8. — Histogram of the fesc,rei limits of our sample. The top panel limits are derived 
assuming a UV-to-LC ratio (/i5oo//70o)ste« = 81 whereas the lower panel show the limits 
derived when taking the UV-to-LC ratio from the best fit model. Limits of fesc,rei < 1-0 
indicate additional attenuation (by either dust or HI) in addition to the attenuation at 1500 
A. Eighteen (of 21) galaxies have limits of fesc,rei < 1-0 and 9 have fesc,rei < 0.5 (when using 
SED fits). 
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Fig. 9. — Best fit SEDs for tliree galaxies. Plotted in eacli panel are the original SED (thin), 
with best-fit reddening applied (thick), and with average IGM absorption applied (blue). Top 
Left: This galaxy (J033237. 07-274617. 3) has no indication of a Balmer break but has a red 
spectral slope so it is fit with a young (30 Myr) and dusty (E(B-V)=0.35) model Top Right: 
A red galaxy (J033244. 16-274729.5) with a prominent Balmer break fit with a very old (3 
Gyr), relatively dusty (E(B-V)=0.25) template. Bottom: A galaxy (J123643.41+621151.6) 
with an SED very similar to typical LBGs (Age = 300 Myr, E(B-V)=0.15). The far-UV 
3a limits are plotted with downward arrows. Note that our ability to detect leaking Lyman 
continuum is largely dependent upon the level of dust attenuation. 
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Fig. 10.— The HST/WFPC2 B-band (£450, left) and far-UV (F150LP, right) images of 
J123652. 69+621355. 3. The orientation and alignment is the same in both images and the 
pointers are in the same location on the sky. The pointers are each 1" in length and point 
to the faint source to the North of the target galaxy. The source is clearly detected in the 
far-UV image. 
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Fig. 11. — The parameter space exc luded by a Monte-Carlo analysis of the combined limits of 
this work with lMalkan et al.l (120031 ). The x-axis is the relative escape fraction and the y-axis 
is the fraction of galaxies which have this escape fraction. The other galaxies are assumed 
to have negligible escape fractions. The shaded regions are excluded(at 3, 2, Icr darker to 
lighter) by the fact that we detect no galaxies in the combined sample of 32 g a laxies . The 
dashed box denotes the approximate parameter space implied by lShapley et al.l (120061 ). The 
horizontal errors to the box are from uncertainties in the determination of the relative escape 
fraction and the vertical errors are Poisson. 



